Dilute Fermi gas in quasi-one-dimensional traps: From weakly interacting fermions 
via hard core bosons to weakly interacting Bose gas 
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We study equilibrium properties of a cold two-component Fermi gas confined in a quasi-one- 
dimensional trap of the transverse size l±. In the dilute limit (nl± <C 1, where n is the ID density) 
the problem is exactly solvable for an arbitrary 3D fermionic scattering length of. When 1±/clf 
goes from — oo to +oo, the system successively passes three regimes: weakly interacting Fermi gas, 
hard core Bose gas and weakly coupled Bose gas. The regimes are separated by two crossovers at 
of ~ ±nij_- In conclusion we discuss experimental implications of these results. 
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Trapped cold atomic gases [1] offer a unique possi- 
bility to create experimentally various many-body sys- 
tems that for a long time have been considered as purely 
theoretical models. In this respect the study of quasi- 
one-dimensional (quasi-lD) systems is especially promis- 
ing since ID many-body problems are frequently ex- 
actly solvable [2]. Currently quasi- ID cold atomic Bose 
gases are indeed extensively studied both experimentally 
[3] and theoretically [4, 5, 6]. The behavior of two- 
component Fermi systems is expected be more diverse 
due to one more "degree of freedom" related to the for- 
mation of two-particle composite bosons (dimers). In 
quasi-3D traps both Bose-Einstein condensation (BEC) 
of dimers [7] and the regime of a crossover from Bardeen- 
Cooper-Schriefer (BCS) superfluidity to molecular BEC 
[8] have been recently reproduced experimentally. While 
theoretical studies of BCS-BEC crossover in 3D (or quasi- 
3D) systems have a long history (see Refs. 9, 10, 11, 12, 
13, 14 and references therein), a similar problem for a 
strongly anisotropic confinement has non been addressed 
up to now. In this paper we present a complete theory of 
"BCS-BEC" transformation from weakly coupled Fermi 
gas to Bose gas of diatomic molecules in quasi- ID traps. 

We consider a dilute two-component Fermi gas con- 
fined in (x, y) plane by a harmonic potential with charac- 
teristic scale lj_ = (mtjjj -1 / 2 . In the present context the 
diluteness means that nl± -C 1, where n is the density 
averaged over transverse directions. As usual, assum- 
ing that the interaction potential V(r) is of short range 
Rw < ii, we employ the standard pseudopotential ap- 
proximation: V(r) = 4 ^° F 5{r) J^r, where ap is the 3D 
scattering length. Parameter l±/ap ranges from — oo to 
+oo, which corresponds an attractive interaction and in- 
cludes the unitarity limit (ap = oo). We will show that in 
the whole range of ap this system is described by exactly 
solvable models. At small negative of {Ix/o-f <C —1) it 
behaves according to Gaudin's solution [15] for an attrac- 
tive ID Fermi gas. At of ~ —nl\ the Gaudin's Fermi gas 
transforms to Girardeau- Tonks gas of hard core bosons 
[16, 17], which, at positive ap > nl\, smoothly matches 
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FIG. 1: Phase diagram of quasi-lD Fermi gas. 
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the regime of Lieb-Liniger Bose gas [18] with weak re- 
pulsion. A schematic phase diagram is shown in Fig. 1. 
Two crossovers in Fig. 1 reflect two different physical phe- 
nomena. The left crossover corresponds to the transfor- 
mation from strongly overlapped BSC-like pairs to well 
defined composite bosons, which is roughly similar to the 
common BCS-BEC crossover in 3D systems. The right 
crossover is inherent for quasi- ID geometry. It is related 
to a change of boson-boson interaction due to a "dimen- 
sional transformation" from quasi- ID to 3D character of 
composite bosons. 

Many-body physics of an attractive Fermi gas is closely 
related to the solution of a two-body problem. The two- 
body scattering problem in an axially symmetric har- 
monic trap has been already addressed by several au- 
thors [4, 5, 19]. Below we briefly discuss properties of 
a two-body bound state that is a localized solution to 
Schrodinger equation for the relative motion. The en- 
ergy of this state should have a form Eq = u>± — Eq, 
where £o > is the binding energy. In the pseudopoten- 
tial approximation the bound state wave function, Xo( r ), 
is expressed in terms of Green's function, G(E, r,r'), of 
a cylindrical harmonic oscillator with the mass m/2: 



X o(r)-AG o (r,0), 



(1) 



Here Go(r, 0) = G{E$, r,0), and A is the normalization 
constant which is given by the expression 
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where £(z,f2) is the zeta function [20]. The binding en- 
ergy Eq is a solution to the following transcendental equa- 
tion 
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(3) 
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Equation (3), which defines a pole in the two-body scat- 
tering amplitude, can be reduced to the form [5] 



ix/o F = -C(l/2,eo/2w±)/v^ 



(4) 



Right hand side in Eq. (4) is a monotonically increas- 
ing function of ft = ea/2cjj_. It diverges as — 1/V20 at 
fl <C 1, crosses zero at O w 0.3, and goes to +oo as v / 2?} 
at £1 3> 1 . This behavior translates to the monotonic de- 
pendence of the binding energy £p on aj?. For small neg- 



ative a F (ii/flF <C —1) we get the result eo 



? /ml\. 



In this regime eo <C which means that the transverse 
motion of particles is confined to the lowest state of spa- 
tial quantization. The bound state is strongly anisotropic 
with the axial size ~ 1 / y/mso ^> l±- In the unitarity 
limit, l±/a F = 0, the bound state wave function becomes 
almost spherically symmetric, while the energy takes a 
universal form Sq = 0.6loj_. In the regime of strong at- 
traction (l±/a F 3> 1) the binding energy approaches the 
usual 3D expression: e ~ l/ma%. 

We also introduce an effective size of the bound state 
ao = l/^/m£oi which has the following behavior in the 
characteristic regions of 1±/clf 

a xsl 2 ± /\a F \ > l±, l±/a F <^-l, (5) 
a ~1.29/_L, lj_/a F = 0, (6) 

ao~flF«ii, l±/a F » 1- (7) 

Let us turn to the many-body problem in the limit 
nlj_ <C 1. We start from the regime of weak attraction: 
a F < and \a F \ <C l± (the region to the left of the shaded 
region in the phase diagram of Fig. 1). In this case the 
transverse motion of each Fermi particle is confined to the 
lowest oscillator's state. Therefore an effective ID many- 
body Hamiltonian is obtained from the original one by 
the simple averaging over transverse directions 



/[ 2 d 2 



(8) 



where i()j(z) are the Fermi field operators, and the ID 
fermionic coupling constant takes the form 



g F = 2\a F \/mi]_. 



(9) 



Solving the ID two-particle problem with the point at- 
tractive interaction of Eq. (9) we find the bound state 
energy eo = m g F /^ = o 2 F /'rrd\, which recovers the cor- 
responding solution to Eq. (4). 

The problem of ID spin-1/2 fermions with an attrac- 
tive interaction, Eq. (8), is exactly solvable [15]. The 



ground state energy E F is given by the set of Gaudin's 
integral equations [15] 
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7^f(k')^, (11) 
(12) 



The solution of Eqs. (10)-(12) is governed by dimension- 
less parameter 7q = mg F /n <~ l/na ~ \a F \/nl\. Since 
7g is a ratio of two small parameters (\a F \/l± and nl± 
respectively) it can take any value in the region of appli- 
cability of Eqs. (8)-(12). 

In the weak coupling limit, 7g <SC 1 (l±/a F -c 
— l/n/j_), we get a weakly interacting two-component ID 
Fermi gas with the energy E F = ir 2 n 3 /24m. To reveal 
the physics behind the strong coupling limit, 7g > 1 
(l_i_/a F <C —l/nl±), it is instructive to rescale all mo- 
menta, k — ► p/2, in Eqs. (10)-(12) and rewrite them as 
follows 
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where f(p) — 2/(p/2), and M — 2m. Obviously, the 
first term in Eq. (13) is the energy of n/2 noninteracting 
bound states (composite bosons). The rest of the set of 
Eqs. (13)-(15) coincides (except for the sign of the inte- 
gral term in Eq. (14)) with a set of Lieb-Liniger equations 
[18] for a gas of interacting Bose particles with the mass 
M and the density n/2. This analogy was already noted 
in the early paper by Gaudin [15]. In the strong coupling 
limit the integral term in Eq. (14) (with a "wrong" sign) 
is irrelevant. Hence in this regime the system behaves 
as a gas of impenetrable composite bosons (Girardeau- 
Tonks gas) . The strong repulsion is a direct consequence 
of Fermi statistics and ID kinematics - the composite 
particles can not go around while the Pauli principle for- 
bids them to penetrate each other. The energy in this 
regime is given by the formula for a one-component Fermi 
gas of the density n/2: Ef = -e n/2 + ir 2 {n/2) 3 /6M '. A 
crossover from the weakly coupled Fermi gas to the Tonks 
gas is located at l±/a F <~ —l/nl±. The description of our 
system in terms of Gaudin's Hamiltonian, Eq. (8), with 
g F of Eq. (9) is valid in the region l±/a F < — C, where 
constant C satisfies the condition l/nl± » C » 1. 

Another asymptotic regime corresponds to a F > and 
Of < (i (the region to the right of the shaded region 
in Fig. 1). This is the regime of small 3D composite 



3 



bosons with the size ao = clf (see Eq. (7)), the binding 
energy Sq = l/map and 3D scattering length as ~ 0.6clf 
[21, 22]. Since nZj_ <C 1, the transverse center-of-mass 
motion of dimcrs corresponds to the lowest harmonic os- 
cillator state. Hence the effective ID bosonic Hamilto- 
nian reduces to the form 
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dz 
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(16) 



with the following ID bosonic coupling constant 

g B = 2a B /ml\ w l.2a F /ml\. (17) 

The energy of n/2 bosons with a repulsive contact in- 
teraction is given by the exact solution due to Lieb and 
Liniger [18] (Lieb-Linigcr integral equations are equiva- 
lent to Eqs. (10)-(12) with the replacement gp — > —gB)- 

A dimensionless parameter that governs the behavior 
of the system in this regime is -jll = MgB/n <~ apjnl\. 
In the region l±/a,F 3> l/nl± (which corresponds to 
111 <C 1) we have a weakly coupled Bose gas with 
the energy E% = -e n/2 + ±g B (n/2) 2 [18]. In the op- 
posite limit lj_/ap <C l/nlj_ ("/ll 3> 1) we again re- 
cover the hard core ID bosons, and thus get the energy 
Eg = -e n/2 + 7r 2 (n/2) 3 /6M. A crossover between 
these two regimes corresponds to lj_/a,F ~ l/nl±. 

Inside the shaded region in Fig. 1 both above asymp- 
totic approaches (Eqs. (8)-(9) and Eqs. (16)-(17)) are not 
applicable. In this region the relative motion two parti- 
cles is neither ID (as assumed in Eqs. (8)-(9)) nor 3D 
(which is the condition for applicability of Eq. (17)). We 
have, however, seen that on either side of the shaded 
region the system behaves a gas of hard core compos- 
ite bosons. Below we prove that this model is also 
valid everywhere inside. The main point is that for all 
l±/a,F S> — l/nl±, which also covers the questionable 
region, the bound state's size a , Eqs. (5)-(7), is much 
smaller than mean interparticle distance 1/n. This al- 
lows us to apply the functional integral approach that 
has been originally developed to describe the molecu- 
lar limit of BCS-BEC crossover in homogeneous systems 
[9, 10, 11]. The generalization of this theory to arbitrary 
quasi-3D trapped systems is given in Ref. 22. Since in 
harmonic traps the center-of-mass and the relative mo- 
tions are decoupled, the general formalism of Ref. 22 is 
directly applicable to the present quasi-lD problem. The 
result of this approach is quite simple and physically ap- 
pealing - to the leading order in a^n <C 1 the system is 
still described by the bosonic Hamiltonian of Eq. (16), 
but gs in general has a more complicated form. In the 
present context it is enough to consider a Born approx- 
imation to the coupling constant gs, which corresponds 
to an "exchange" process shown in Fig. 2 [13, 22]. An 
analytic expression for this diagram takes a form 

f 4 

9B=K J t B (Pl, P2, P3, Pi) J\_^0fi{Pi)d 2 Pi, ( 18 ) 




FIG. 2: Diagram for the Born scattering amplitude of two 
composite bosons. 



where p = (x, y) is the transverse coordinate, 
$o,o(p) = l^ 1 yj2/TTe~( p / [± -^ 2 is the wave function for 
the transverse motion of the center-of-mass, and Ao = 
JV(r)xo{r)d 3 r = A (see Eqs. (1), (2)) is the boson- 
fermion "vertex" [22]. The four-point function is is de- 
fined as follows 

t B = f 77^Gu,k(PUP2)GZ,k{P2>P3) 



, (2tt)2- 

X0u;,k(P3)P4)0w,fc(P4,Pl), 



(19) 



Here Gu,k{p, p') is the one particle Green's function for 
a harmonically confined Fermi gas with the chemical po- 
tential fj, = oj± — £o/2 
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(20) 



where 4> n ,m{p) are the eigen functions of a 2D harmonic 
oscillator with frequency uj±. Reducing all integrals in 
Eqs. (18), (19) to a dimensionless form we find that gs 
has the following general structure 



9 b = 9{eo/uA_)/mh 



(21) 



Dimensionless function g(x) increases as 6\/x at x <C 1, 
reaches a maximum at i ~ 1, and goes to zero as 
4/^ at x ^> 1. Using the known dependence of e on 
Ij_/cif we get the required function gs{aF) in all char- 
acteristic regimes: gB ~ \aF\/ml\ if Zj_/|of| 3> 1, and 
gs ~ l/ml±, if l±/ 1 Of | < 1 [23]. Renormalization of 
the Born scattering amplitude does not alter this gen- 
eral behavior since it can change only numerical coeffi- 
cients in the above formulas. Thus, everywhere in the 
region ?j_/|(If| <C l/nl± the strong coupling condition 
mgB/n S> 1 is fulfilled. This proves the correctness of 
the complete phase diagram in Fig. 1 . 

The above formal results lead to the following physical 
picture of the quasi- ID "BCS-BEC" transformation that 
consists of two well defined steps. The first step is the 
"Fermi-Bose" crossover at l±/a F <~ — l/nl±. On the left 
of this point we have a "BCS"-like state with strongly 
overlaped pairs, while everywhere on the right the sys- 
tem is composed of well defined composite bosons. The 
further physical changes in the system are related to the 
evolution of an effective boson-boson interaction. In the 
region —l/nl± < l±/ap < C (from the left crossover to 
the right end of the shaded region in Fig. 1) the bosons 
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are of quasi- ID form. Their transverse size equals to the 
transverse size l± of the confining potential. Therefore 
two colliding bosons can not avoid each other while their 
mutual penetration is forbidden by Pauli principle. As 
a result we have Girardeau- Tonks regime. In the region 
C > l±/a,F > l/nl± the composite bosons take the com- 
mon 3D form with radius cio ~ aj?. However the trans- 
verse confinement of their center-of-mass motion is still 
strong - the system remains in the hard core gas regime. 
It should be outlined that the microscopic reason for the 
Tonks gas regime is quite different on opposite sides of 
the shaded region. Finally when the size of bosons be- 
comes small enough (do ~ ap < nlj_) the system enters 
the regime of weakly interacting Bose gas. It is worse 
mentioning that in contrast to 3D case nothing special 
happens in the unitarity limit = oo that is simply 
the middle of the Giradeau- Tonks regime. In fact this 
regime is a quasi-lD realization of a universal behavior 
when all thermodynamic characteristics are independent 
of the fermionic scattering length. An important fea- 
ture of quasi- ID systems is that the universal behavior 
extends to a wide region around the point ay = oo. 

Experimentally quasi-lD atomic gases are normally 
produced in highly elongated needle-shaped traps [3] 
while the scattering length clf can be tuned using the 
Feshbach resonance technique [7, 8] . In such systems the 
double-crossover structure of the phase diagram can be 
observed via changes of the density profile n(z). If the 
axial confinement (with a scale l z 3> l±) is semiclassical 
we can calculate n(z) using the local density approxi- 
mation. In the regime of weakly interacting Gaudin's 
gas we get the result uq(z) = (2/nl z )y/N — z 2 /l 2 , which 
corresponds to noninteracting fermions (here N is the to- 
tal number of particles). The density distribution in the 
Girardeau- Tonks regime is given by the similar expres- 
sion 71GT 

Hence when we scan through the first crossover at 
lj_/aF ~ — l/n(0)Zj_ <~ — l the cloud shrinks down 
while preserving the Fermi-like shape. The density dis- 
tribution remains unchanged until the next crossover at 
1±/cif ~ 7~^//v * s reacne d- This is the universal or uni- 
tarity limited regime as it includes the point of = oo. 
In the region Ij_/clf > j the system gradually trans- 
forms to a weakly coupled Bose gas with a parabolic den- 
sity profile tib(z) = ras(0)(l — z 2 /Z<^ F ), where n#(0) = 
a^il^ZTp/l 2 ) 2 and Z TF = (3a B ^A r /4?i) 1/3 . The evo- 
lution of n(z) can be visualized in terms of the follow- 
ing two-step dependence of the mean-square size, (z 2 ) = 
N^ 1 J z 2 n(z)dz, on lj_/a F 



{4 z ' a F lj_VN 

& N > ^ \^<TZ7W (22) 

1/2 1 3ha B N \ 2 ^ 3 l± . l z 

5 l 2 V il l J ' «f IiViV' 

The most characteristic feature of quasi- ID Fermi 



gases is a long, well pronounced plateau with (z 2 ) w 
\l 2 z N in the region Z_L/|a F | < j-^- (see Eq. (22)). This 
should be contrasted to the plain single-crossover struc- 
ture of (z 2 ) in 3D traps, which has been observed in 
Rcf. [8]. Therefore the observation of such a plateau 
around the Feshbach resonance in a highly elongated 
trapped Fermi system unambiguously indicates a real- 
ization of the quasi-lD regime. 
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